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ABSTRACT 

We investigate stability of a gravitationally collapsing iron core against non-spherical 
perturbation. The gravitationally collapsing iron core is approximated by a similarity 
solution for dynamically collapsing polytropic gas sphere. We find that the similarity 
solution is unstable against non-spherical perturbations. The perturbation grows in 
proportion to (t — to) _<J while the the central density increases in proportion to (t — 
to) -2 - The growth rate is a = 1/3 + £(7 — 4/3), where 7 and £ denote the polytropic 
index and the parameter £ of the spherical harmonics, Y™(0, tp), respectively. The 
growing perturbation is dominated by vortex motion. Thus it excites global convection 
during the collapse and may contribute to material mixing in a type II supernova. 

Subject headings: gravitation — hydrodynamics — instabilities — stars: supernovae 

1. INTRODUCTION 

Supernova remnants are appreciably aspherical and globally asymmetric. The asymmetry of 
ejecta indicates that the supernova explosion is highly non-spherical and contains non-radial flow 
(see, e.g., the review by Goldreich, Lai & Sahrling 1997 and references therein). The non-radial 
flow in supernova explosion is suggested also from x-ray and 7-ray observations of SN1987A; one 
cannot explain early detection of x-rays and 7-rays from SN1987A without invoking large scale 
mixing (see, e.g., the review by Bethe 1990 and the references therein). The origin of high velocity 
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pulsars may also be ascribed to asymmetry of supernova explosion (see, e.g., Burrows & Hayes 1996 
and references therein). 

Asymmetry is amplified by the Rayleigh- Taylor instability in the early phase of supernova 
explosion (see, e.g., Falk & Arnett 1973). According to detailed numerical simulations, the growth 
of the Rayleigh- Taylor instability accounts for matter mixing inferred from observations of SN1987A 
if there exists an appropriate seed of the asymmetry (Arnett, Fryxell, & Miiller 1989; Hachisu et 
al. 1990; Miiller, Fryxell, & Arnett 1991; Fryxell, Arnett, & Miiller 1991; Nagataki, Shimizu, k 
Kato 1998). Thus it is worth to consider the origin for the seed of asymmetry. 

Goldreich et al. (1997) discussed possible instabilities during the core collapse as a source of 
the seed. They suggested possibility that dynamically collapsing core might be unstable against 
non-radial perturbation and proposed to study the stability of the similarity solution of Yahil 
(1983). Applying the polytropic equation of state, P = Kp" ( , to pre-supernova core, he obtained 
a similarity solution describing collapse of a spherical iron core. In this paper we show that his 
similarity is indeed unstable against a vortex mode. The method of stability analysis is essentially 
the same as that of Hanawa & Matsumoto (2000) who investigated the bar mode instability during 
collapse prior to protostar formation. Velocity perturbation dominates over density perturbation 
in the vortex mode while both of them have similar amplitude in the bar mode. The growth rate of 
the vortex mode depends on the wavenumber, £, of the spherical harmonics, Y™(0, (p). The vortex 
mode grows in proportion to \t — to\~ a , where to and a = 1/3 + £ (7 — 4/3) denote the epoch of 
protoneutron star and growth rate, respectively. 

In §2 we review the similarity solution of Yahil (1983) for our stability analysis given in §3. 
We discuss the mechanism of the vortex mode in §4 and implications to type II supernova in §5. 
We show the asymptotic behavior of the similarity solution and perturbation in the region very far 
from the center in Appendix. 



2. SIMILARITY SOLUTION 



For simplicity we consider gas of which equation of state is expressed by polytrope, 

P = Kp\ (1) 

where P and p denote the pressure and density, respectively. The hydrodynamical equations are 
then expressed as 

^ + V-(pv) = (2) 

and 

d 

g^(pv) + VP + V- (pv (g>v) + pV$ = 0, (3) 
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where v and denote the velocity and gravitational potential, respectively. The gravitational 
potential is related with the density distribution by the Poisson equation, 

A$ = 4vrGp, (4) 

where G denotes the gravitational constant. 

For later convenience, we introduce the zooming coordinates of Bouquet et al. (1985) to 
solve equations (1) through (4). The zooming coordinates, (£, r), are related with the ordinary 
coordinates, (r, t), by 

( *) = [ c o\t-t \ ] , (5) 
V r / V - In |1 - t/t \ J 

where cq denotes a standard sound speed and is a function of time t. The symbol, to, denotes an 
epoch at the instant of the protoneutron star formation. The density in the zooming coordinates, 
g, is related with that in the ordinary coordinates, p, by 

q(x, t) = ^G P {t - t ) 2 . (6) 

We define the standard sound speed, Co, so that it denotes the sound speed at a given t when 
g = 1. Thus it is expressed as 

n(1-7)/2u _ , ,1-7 



Co = y'yK (Ai^Gy- ' ' \t — to\ ■ (7) 
The pressure in the zooming coordinates, p, is related with the that in the ordinary coordinates, 

AttG 

Substituting equations (6) and (8) into equation (1), we obtain the polytrope relation in the zooming 
coordinates, 

P = - ■ (9) 



P = —P(t ~ to) 2 • (8) 



7 

The velocity in the zooming coordinates, u, is defined as 

u = y + (2 - 7 ) / . (10) 
Co Co |t — *o I 

This velocity denotes that with respect to the zooming coordinates, and includes the apparent 
motion, the last term in equation (10). The gravitational potential in the zooming coordinates, (p, 
is related with that in the ordinary coordinates, <E> by 

In the zooming coordinates, the hydrodynamical equations are expressed as 

|f + V r (qu) = (4 - 37) Q , (12) 
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d 

— (gu) + V f • {gu ® «) + V 5 p + ^V e = (2 - 7) (7 - 1) ^ + (7 - 5 7 ) £U , (13) 



and 



= q (14) 

for t < to. The symbols, and A^, denote the gradient and Laplacian in the £-space, respectively. 

Assuming stationarity in the zooming coordinates [d/dr = 0) and spherical symmetry 
(d/d6 = d/dip = 0), we seek a similarity solution. Under these assumptions equations (12), 
(13), and (14) reduce to 

du r u r dg 2u r 

-^7- + — = (4-37 — , 

9$, g d£ £ 



and 



du r 1 (dp\ dg d(f> 1U , /Q . 

®4> 1 / ,1 QUr 



(15) 
(16) 



where £ = |£|. After some algebra we can rewrite equations (15) and (16) into 



,7-1 



V) ( I ) = e 



gu r 



4-37 



+ (2 - 7 ) (7 - IK + (7 - IK + 



2u r 2 ' 



(17) 



(18) 



and 



,7-1 



u r 2 ) 



gu r 



4-37 

+ (4-3 7 )^- 1 



(2 - 7) (7 - l)£u r - (3 - 2 7 ) u 2 
,7-1 



2« r 



(19) 



Equations (18) and (19) are singular at the sonic point, u r 2 = £> 7_1 . We obtain the the similarity 
solution by integrating equations (18) and (19) with the Runge-Kutta method. In the numerical 
integration we used the auxiliary variable of Whitworth Sz Summers (1985), s, defined by 



— = g 1 - U r . 

as 

Using equation (20), we rewrite equations (18) and (19) into 



(20) 



dg 
ds 



gu r 



4-37 



+ (2 - 7X7 - ik + (7 - IK + 



2u r 



and 



(21) 



du r 
ds 



gu r 



4-37 
+ (4 - 3 7 ) (f- 1 



(2 - 7) (7 - 1) tu r - (3 - 2 7 ) u r 



2u 



V 1 



(22) 
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respectively. 

Similarity solutions exist for 7 < 4/3. Figure 1 shows the similarity solution for 7 = 1.3. The 
solid curves denote g while the dashed curves denote the infall velocity, — v r = — u r + (2 — 7) £. 
These solutions are the same as those obtained by Yahil (1983) and Suto & Silk (1988). They have 
the asymptotic forms of 



Qc - 



2-7 

Qc 



(23) 



and 



u 



r 7 ) ? + 3 



The value of g c is 22.04 for 7 = 1.3. 



j -7) e 5 + o(e) 



H ■ 



(24) 



3. VORTEX MODE 

In this section we consider a non-spherical perturbation around the similarity solution. The 
density perturbation is assumed to be proportional to the spherical harmonics, Y™{6, (p). Then 



the density and velocity are expressed as 

Q = Qo + 5Q(Oe° T Y e m (9,<p), (25) 

u r = u r0 + 6u r (Oe aT Y e m (0,<p), (26) 

ue = Su e (0-^j^Yr(e^), (27) 

u v = MO (<+ ^^ (^). (28) 

4> = 0o + 64>(Oe aT YP(9,<p) , (29) 

where the symbols with suffix, 0, denote the values in the similarity solution and the symbols 
with the symbol, 5, denote the perturbations. Substituting equations (26) throughout (29) into 
equations (12), (13), and (14), we obtain the perturbation equations, 

( ff + 37 - 4)Sg + lA[e 2 (^o<5n r + u r0 Sg)} - I ^ = 0, (30) 

(a + 2 7 - 3) 5u r + ^ (u r0 6u r ) + ^ (j^j + 6F = ' (31) 

(a + 2 7 - 3)5u e + ^ ± + l -±± (-^ + 8<ty = , (32) 

— (50 = <5r , (33) 
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and 



d _ 2 5r £{£ + 1) ^ . 
^ST = -— + -^-Ls<I, + Sq. 



(34) 



These perturbation equations have singularities at the origin (£ = 0), the sonic point [(u r o) 2 — 
1 ], and the infinity (£ = +oo). These perturbation equations are the same as those of Hanawa 
& Matsumoto (2000). 

To obtain the boundary condition at the origin we use the Taylor expansion of the perturbation. 
In the following we use the notation, 



Sg 



go 



2-7 

Su r 
8uq 



and 



h = \ e + Me +2 + o(^ +4 ). 



tt + i\ 



Substituting equation (38) into equation (33) we obtain 

ST = l\o£ e - 1 + (e + 2)A 2 £' + 1 + 0(£ e + 3 ) . 



(35) 

(36) 
(37) 

(38) 
(39) 



Note that not only the leading terms but the second lowest terms are taken into account in equations 
(35) through (39). 

Substituting equations (35) through (39) into equations (30) through (34) we derive conditions 
for ao, tt2, Po, P2, £q, £2, Aq, and A2. From equation (30) we obtain 



it + 1) «o 



'■Po 







(40) 



Equation (40) ensures that the terms proportional to £^ -2 vanish in equation (30). The terms 
proportional to £ £ ~ 2 vanish in equation (34) at any condition. From equations (31) and (32) we 
obtain 

'4 

- 7 



a + 2 7 - 3 + I 



and 



a + 2 7 -3 + ^(--7 



"o + (■ (£0 + A ) = , 



Po + {I + 1) (£ + A ) = 



(41) 



(42) 



respectively. These equations ensure that the terms proportional to ^ 1 vanish in equations (27) 
and (28). Similarly we obtain 



2- 7 



7 



£0 + Qc[{t + 3)q 2 - efo] 

0, 



Qc ~ 3 ) [(* + 3)«o - £0o 



(43) 
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and 



[(£ + 2) + 3) - £ {£ + 1)] A 2 



2-7 



from the condition that the terms proportional to £ £ vanish and 



a + 2 7 - 3 + {£ + 2) 



- 7 



a 2 + + 2) 



1-7 

c 

15 



Qc ~ 



(44) 



- 7 a 



+ {£ + 2) (e 2 + A 2 ) = , 



(45) 



and 



a + 2 7 - 3 + {£ + 2) 



7 



ft + 



1-7 

c 

Is - 



7 A) 



+ (£ + 1) (e 2 + A 2 ) = 

from the condition that the terms proportional to £ +1 vanish. 
From equations (40), (45) and (46) we obtain 



(46) 



7 



[(e + i) a 2 - (e + 2) &] = o . 



(47) 



(48) 



a + 2 7 - 3 + (£ + 2) 

This condition is equivalent to the condition that either of 

(I + 1) q 2 - (I + 2) 2 = 
and 

,_| + /( 7 -|) (49) 

is fulfilled. The bar mode found by Hanawa Sz Matsumoto (2000) fulfills equation (48). We seek 
the other mode that fulfills equation (49). In the following we call the latter the vortex mode. 

Since the perturbation is linear, the solution can be arbitrarily scaled. To normalize the 
solution we take qq = 1 in this paper. Then we obtain 



0o 



Ao 



£o = -ZqI- 1 [{Z + 3)a 2 - £ 2 ] + ( q, 



i + 1 

t 

2 /4 



Ur 7 



2-7 

Qc 



2 (2£ + 3) 



and 



e 2 



1-7 

Qc 

15 



7 



\2 • 



(50) 
(51) 
(52) 
(53) 

(54) 
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We obtained an eigenfunction of this vortex mode numerically by the following procedures. 
First we integrated equations (21) and (22) to obtain the similarity solution for a given 7. Second 
we obtained three linearly independent solutions for the perturbation around the origin by using 
the Taylor series expansion, equations (35) through (39). Third we integrated the three linearly 
independent solutions from the origin toward the sonic point numerically with the Runge-Kutta 
method. Forth we obtained two linearly independent solutions satisfying the boundary conditions 
both at the origin and sonic point by taking linear combinations of the numerically integrated 
solutions. Finally we integrated the two linearly independent solutions from the sonic point toward 
the infinity and obtained an eigenfunction satisfying all the boundary conditions as a linear combi- 
nation of them. The eigenfunction should have an infinitesimal small amplitude at the infinity as 
shown in Appendix. 

Figures 2, 3, and 4 denote the numerically obtained eigenfunctions of I = 1, 2, and 3, respec- 
tively. The polytropic index is set to be 7 = 1.3 in the figures. The eigenfunctions are normalized 
so that the radial velocity perturbation is 5u r = + O (^ + 1 ) near the origin. The non-radial 

velocity perturbation, Su$, changes its sign around £ « 1. The density perturbation is small. It 
should be proportional to (7 — 4/3) [see eq. (51)]. 

The vortex mode of i = 1 is different from the ghost mode of i = 1 (Hanawa & Matsumoto 
1999, 2000). The vortex mode denotes circulation streaming back through the core surface while 
the ghost mode denotes misfit of the gravity center to the coordinate center (Hanawa & Matsumoto 
1999). The ghost mode has no vortex (V^X(5w = 0). The vortex mode of I = 1 has the growth 
rate of a = 7 — 1 while the ghost mode has that of a = 2 — 7. 

Figure 5 shows the density and velocity perturbation in the cross section. The contours denote 
the iso-density curves of g = 10, 3, 1, 0.3, and 0.1 while the central density is g c = 22.04. The 
arrows denote the velocity vectors. We obtained the density and velocity by adding the eigenmode 
of (£, m) = (2, 0) on the similarity solution for 7 = 1.3. The infall is a little faster in the x-direction 
than in the z-direction. Figure 6 is the same as Figure 5 but only the velocity perturbation is shown 
by the arrows in Figure 6. It shows that this eigenmode is an vortex flow in the meridional plane. 



In this section we discuss the mechanism for the growth of the vortex mode. We consider the 
conservation of vorticity to elucidate the growth mechanism. Taking rotation of the equation of 
motion, we obtain, 



4. 



GROWTH MECHANISM 




(55) 



where 



ft = VXD , 



(56) 
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since the rotation of the pressure force and that of gravity vanish. Similarly the conservation of 
the vorticity is expressed as 



where 

Since V^X uq = 0, we obtain 

d_ 



V ? x (ux«) + (3 - 27)0; , 
uj = V^xti. 

5u = V 5 x (tt x 5u) + (3 - 27) 5u , 



(57) 
(58) 

(59) 



for the perturbation of the vortex, 5uj. Substituting the equation (24) into equation (59) we obtain 



dr 



5u 



7 



a 



3 J <91n£ 



+ 



8uj , 



(60) 



near £ = 0. We can derive equation (49) from equation (60) since 5u oc ^ near £ = in the 
vortex mode. We can derive also the growth rate of the spin-mode, a = 1/3, from equation (60) 
since 5u> = £° near £ = in the spin-up mode. The growth of the vortex mode as well as that of 
the spin-up mode is subject to the conservation of the vorticity. 



5. IMPLICATION TO TYPE II SUPERNOVA 

As shown in the previous sections the vortex mode grows in proportion to \t — to|~ CT - In other 
words it grows in proportion to p^ 2 since p c oc (t — to)~ 2 . This growth rate is for the growth of the 
relative amplitude, i.e., that for Su = 5v/c$. The growth of the anisotropic velocity is proportional 
to 

5v oc c Pc /2 oc pi , (61) 

where 

✓ -i + £±i(,-i). 

Suppose that an iron core during implosion can be well approximated by a polytrope of 7 = 1.3. 
Then the velocity perturbation of I = 1 grows by a factor of 30 while the density increases from 
Pc = 10 9 g cm 3 to 10 14 g cm 3 . Similarly that of I = 2 grows by a factor of 25 during the same 
period. 

If the I = 1 mode is amplified, the central core, i.e., the protoneutron star has a bulk velocity 
relative to the envelope. This might explain a run-away pulsar from its natal nebula. If the 
I > 2 mode is amplified, the anisotropic velocity will cause global mixing during the supernova 
explosion (t > to). This may explain heavy element exposure earlier than expected from a spherical 
symmetric model. 
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It should be also noted that the vortex mode of (£, m) = (2, 0) has velocity field similar to 
that of the Eddington-Sweet meridional circulation. If a pre-supernova star rotates slowly, the 
meridional circulation can be a seed of the vortex mode amplified during the implosion phase. 
Convection during Si-burning may also be the seed for the vortex mode. 

We thank Shigehiro Nagataki, Katsuhiko Sato, and Shoichi Yamada for helpful discussion. This 
research is financially supported in part by the Grant-in- Aid for Scientific Research on Priority Areas 
of the Ministry of Education, Science, Sports and Culture of Japan (No. 10147105, 11134209). 



A. Asymptotic Behavior around the Infinity 



In this appendix we derive asymptotic forms of perturbations around the similarity solution 
for a collapsing gas sphere. In the region of £ S> 1 the similarity solution has the asymptotic form 
of 

,ocr 2/{2 - 7) , (ai) 

and 

K - (2 - 7 )£] oc £(i-7)/(2- 7 ) . ( A2 ) 

See Yahil (1983) and Suto & Silk (1988) for the derivation. 

As a boundary condition we assume that the relative density perturbation, 5g/ go, is vanishingly 
small at infinity, £ = oo. After some algebra we obtain the asymptotic relations, 

Sq 



Q0 



oc £ 



-<r/(2- 7 ) 



5u r (X r (- + 7-l)/(2-7) ) 

5u e oc r (-+7-i)/(2-7) 5 

OC r (a-27 + 2)/(2-7) ; 



and 



{a - 2 7 + 2)(<7 - 3 7 + 4) 



- 1(1 + 1) 



r 2 Sg . 



(A3) 

(A4) 
(A5) 
(A6) 

(A7) 



(2 - 7 )2 

When we derive the above relations, we use equations (Al) and (A2). See also Hanawa & Mat- 
sumoto (2000) for the derivation. 
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Fig. 1. — The similarity solution is shown for a collapsing polytropic gas sphere of 7 = 1.30. The 
solid curves denote the density, f?(£), and the infall velocity, —v r (£) = (2 — 7) £ — u r (£). The 
latter is normalized by the sound speed at the center. 




Fig. 2. — The eigenfunction of £ = 1 mode is shown as a function of £. The solid curve denotes the 
radial velocity perturbation, 5u r . The dashed and dotted curves denote 5ug and dp, respectively. 
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Fig. 3. — The same as Figure 1 but for t = 2 mode. 




Fig. 4.— 



The same as Figure 2 but for i = 



3 mode. 
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Fig. 5. — This cross section shows a dynamically collapsing iron core suffering the vortex mode of 
1 = 2. The contours denote the isodensity curves of g = 10.0, 3.0, 1.0, 0.3, and 0.1. The arrows 
denote the velocity, v, in the x — z plane. 




x 



Fig. 6. — The same as Figure 4 but only the velocity perturbation is shown by the arrows. 
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